function x = sitcol(A, b)
% rezolvare orientata pe coloane a sistemului inferior triunghiular
% Intrari: 
%		A=matrice sistem
%	        b=vector termeni liberi
% Iesiri: 
%		  x=solutia sistemului

	[n,n] = size(A);
	x = b;
	for k = 1 : n
		x( k ) = x( k ) / A(k, k);
		x( k+1 : n ) = x( k+1 : n ) - x( k ) * A( k+1 : n, k );
	endfor

endfunction